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The two-dimensional cluster state, a universal resource for measurement-based quantum compu- 
tation, is also the gapped ground state of a short-ranged Hamiltonian. Here, we examine the effect 
of perturbations to this Hamiltonian. We prove that, provided the perturbation is sufficiently small 
and respects a certain symmetry, the perturbed ground state remains a universal resource. We do 
this by characterizing the operation of an adaptive measurement protocol throughout a suitable 
symmetry-protected quantum phase, relying on generic properties of the phase rather than any 
analytic control over the ground state. 



I. INTRODUCTION 

A quantum computer relies on quantum entanglement 
to achieve computational speedups. In the traditional, 
circuit-based model for quantum computation, the re- 
quired entanglement is built up throughout the course of 
the computation through application of entangling gates 
coupling two or more qubits at a time. Alternatively, 
in the model of measurement-based quantum computa- 
tion (MBQC) pfl[2]) universal quantum computation is 
achieved solely through single-particle operations (specif- 
ically, single-particle measurements) on a fixed entangled 
resource state, independent of the quantum algorithm be- 
ing performed. 

Since the initial discovery that the 2-D cluster state 
is a universal resource for MBQC [T], much effort has 
been devoted to characterizing other universal resource 
states. Many of the universal resource states so far iden- 
tified [TJ |3HB] have been projected entangled pair states 
(PEPS) [7] of small bond dimension. The tensor net- 
work structure of these states facilitates the analysis of 
measurements, which might otherwise be an intractable 
problem. Another advantage of such states is that un- 
der appropriate conditions [5], they are unique (possibly 
gapped) ground states of local frustration-free Hamilto- 
nians on spin lattices. This suggests a method of con- 
structing the resource state by cooling an appropriate 
interacting spin system ^] . 

However, if we wish to adopt this viewpoint of the 
resource state for MBQC as the ground state of a quan- 
tum spin system, it would be too restrictive to confine 
ourselves to states in which the effect of measurements 
can be determined analytically from the tensor-network 
structure. A generic local Hamiltonian, or even an arbi- 
trarily small generic local perturbation to a PEPS parent 
Hamiltonian, will not have such a property. Therefore, 
it is desirable to develop an understanding of MBQC in 
ground states of spin systems that does not rely on ana- 
lytic control of the ground state. For this reason, there 
has been an interest in relating MBQC to forms of quan- 
tum order which, as parameters of the Hamiltonian are 
varied, can disappear only at a quantum phase transition 

niHsi- 

In this paper, we will use such a connection between 



MBQC and quantum order to give a precise characteri- 
zation of the operation of MBQC in the ground states of 
a large class of perturbations to the 2-D cluster model. 
This will allow us to give a rigorous proof that such 
perturbed ground states remain universal resources for 
MBQC provided that the perturbation is sufficiently 
small. Our proof relies in part on an extension of the 
the relationship introduced in [T3] between MBQC and 
symmetry-protected topological (SPT) order |14H16j . a 
form of quantum order characterizing quantum systems 
which cannot be smoothly deformed into a product state 
while a certain symmetry is enforced. If the perturbation 
to the 2-D cluster model respects an appropriate sym- 
metry, then the perturbed ground state will still possess 
non-trivial SPT order, and we will show that this gives 
us sufficient information about the ground state to char- 
acterize the implications of the perturbation for MBQC. 
Our result therefore holds independently of any analytic 
solution for the perturbed ground state. 

Our proof of universality is in the same spirit as |17) . 
There, it was shown that, whereas measurements on the 
cluster state simulate quantum circuits, measurements on 
a noisy cluster state simulate the same circuits, but with 
added noise. Here, our task is complicated by the highly 
correlated nature of the "errors" in the resource state 
that from result a change in the Hamiltonian. Neverthe- 
less, we will show how to exploit the additional structure 
resulting from SPT order to establish an effective noise 
model for ground states of appropriate perturbed cluster 
models. Therefore, universal quantum computation can 
be achieved (for sufficiently small perturbations, corre- 
sponding to sufficiently weak noise in the effective circuit 
model) by choosing a measurement protocol which sim- 
ulates a fault-tolerant quantum circuit. The universality 
is then a consequence of the threshold theorem |18j for 
fault-tolerant quantum computation with noisy quantum 
circuits. 



A. Summary of results 

Our ultimate goal in this paper is to prove the uni- 
versality for a MBQC of a class of perturbations of the 
2-D cluster state. However, in order to reach this goal, 
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most of this paper will be devoted to a further elucida- 
tion of the relationship between SPT order and MBQC. 
For simplicity of presentation, we will first explore this 
relationship in one-dimensional systems. It has already 
been shown that in a class of quantum phases charac- 
terized by SPT order, the structure implied by SPT or- 
der leads to the perfect operation of the identity gate in 
MBQC [13] ■ Here, we consider the 1-D cluster model, 
which lies in the simplest of the SPT phases considered 
in |13j . and characterize the operation of non-trivial (i.e., 
not the identity) gates in the presence of a perturbation 
which respects the symmetry protecting this SPT phase. 
We obtain the following: 

Theorem 1 (Effective noise model in one dimension). 

Consider a measurement protocol which in the exact 1-D 
cluster model would simulate a sequence of gates. In the 
perturbed resource state, the same measurement proto- 
col simulates the same gate sequence, but with additional 
noise associated with each non-trivial gate. So long as 
the non-trivial gates are sufficiently separated from each 
other by identity gates, this effective noise has no corre- 
lations between different time steps, i.e. it is Markovian. 

The proof of Theorem[T]will be divided into two stages. 
First, in Section|n]we will establish Theorem[l]for ground 
states which are pure finitely- correlated states (pFCS), 
a special case of matrix-product states (MPS). For such 
states, both the manifestations of SPT order [T51 US], 
and the effect of measurements [3] can be understood 
straightforwardly in terms of the tensor-network struc- 
ture. The ideas leading to Theorem [T] can thus be under- 
stood most directly in this context. Second, in Section 



III we will prove Theorem [T] for arbitrary ground states 
within the SPT phase. 

The extension of these ideas to the 2-D cluster model 
will be considered in Section \FV\ We will construct an ap- 
propriate symmetry group, such that the following result 
is satisfied for symmetry-respecting perturbations. 

Theorem 2 (Effective noise model in two dimensions). 
Consider a measurement protocol which in the exact 2-D 
cluster model would simulate a sequence of gates. In the 
perturbed resource state, the same measurement proto- 
col simulates the same gate sequence, but with additional 
noise associated with each gate. So long as the non-trivial 
gates are sufficiently separated from each other by identity 
gates, this effective noise has no correlations between dif- 
ferent time steps, or between different gates taking place 
at the same time step, i.e., it is local and Markovian. 

Combined with the existing results on fault tolerance 
in the circuit model [18] , Theorem [2] will imply the main 
result of this paper: 

Theorem 3. For sufficiently small symmetry-respecting 
perturbations, the perturbed ground state remains a uni- 
versal resource for measurement-based quantum compu- 
tation. 



II. THE EFFECTIVE NOISE MODEL 
CONSTRUCTION: FINITELY-CORRELATED 
STATES 

In this section, we will prove our effective noise model 
result, Theorem [TJ for a restricted class of ground states 
of infinite one-dimensional chains. Specifically, we con- 
sider pure finitely- correlated states (pFCS) |191 120j . A 
pFCS can be considered as the thermodynamic limit 
of the translationally-invariant matrix-product states 
(MPS) |*jv) generated by the MPS tensor A, on finite 
chains of N sites with periodic boundary conditions, e.g. 
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(here and throughout this paper, we use a graphical no- 
tation to represent the contraction of tensors, e.g. see 
[Him]). The MPS tensor A must satisfy an additional 
condition known as injectivity, which is related to the ex- 
ponential decay of the correlation functions; each of the 
finite-chain states |^jv) (for sufficiently large N) is then 
the unique gapped ground state of a local frustration-free 
Hamiltonian [2"3] , 

We have several motivations for considering this class 
of ground states. First, it is widely believed that pFCS 
capture the essential physics of gapped ground states 
of infinite one-dimensional translationally-invariant spin 
chains in general. (Note that, although the theorem re- 
garding the efficient approximation of ground states of fi- 
nite spin chains by MPS [2H[5S] could be regarded as sup- 
porting this belief, we cannot use this theorem to draw 
any rigorous conclusions for our purposes here, since it 
does not hold that the MPS tensor A can be kept fixed 
independently of the system size for a constant accuracy.) 
Second, the ideas leading to our effective noise model re- 
sult find their simplest and most physically meaningful 
expression in this context. Finally, the proof presented 
here will play a dual role in our paper, as it can also be 
applied to arbitrary quantum states, provided that they 
satisfy a few extra criteria in common with pFCS. Thus, 
in order to establish the effective noise model result for 



general ground states, which we do in Section III it will 
suffice to provide a separate proof of these criteria. 

The outline of this section is as follows. We begin in 
Sees. II A and II B by reviewing the properties of the 1-D 
cluster model and the nature of the SPT phase in which 
it is contained. In Sec. II C we review the results of [T3] 
regarding the structure shared by pFCS ground states 



throughout the whole SPT phase. In Sec. II D we prove 
a key result: the standard adaptive measurement proto- 
col acting on a ground state in the phase is equivalent to a 
non-adaptive dual process acting on a 'topologically dis- 
entangled' version of the ground state, which we refer to 
as the dual state. In Sec. |IIE| we give a characterization 
of the dual process in the case that the original resource 
state is the exact cluster state. Finally, in Sec. II F we 
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U(x) X X X X 

U(z) = X X X X 

FIG. 1. The generators of the on-site x Z2 symmetry of 
the 1-D cluster model. Each dot is a qubit, and the shaded 
areas constitute two-qubit sites. 



exploit the short range of the correlations in pFCS to 
construct the effective noise model for any pFCS ground 
states within the SPT phase, establishing Theorem [l] for 
the case of pFCS ground states. 



A. The 1-D cluster model in the absence of 
perturbations 

Here we recall the properties of the 1-D cluster model 
in the absence of perturbations. The Hamiltonian is 

H = — Zj-iXjZi+i, (2) 



as a linear operator (acting on states from the right), 
where 

0- 

denotes the tensor obtained by complex conjugation from 
the rank-1 tensor corresponding to the state 

The MPS representation for the cluster state plays a 
crucial role in the correlation space picture [3J for the op- 
eration of the cluster state as a quantum computational 
wire [27j . When a projective measurement is performed 
on a site, giving the outcome \ip), this is interpreted as 
inducing an evolution A[ip] on a "correlation system". 
In the case of the cluster state, for any qubit rotation U 
about the x- or z-axis, one can find a product basis {|a)} 
for a two-qubit site, such that 

A c [a] = B a U, (4) 

where B a is an outcome-dependent unitary byproduct 
operator. This byproduct can be accounted for by ad- 
justing the basis for future measurements depending on 
the outcome of the current one. 



where Xi denotes the Pauli X operator acting on the i- 
th site, and similarly for Z{ . With appropriate boundary 
conditions, the system has a unique ground state (the 
cluster state), and an energy gap of 2, independent of 
the system size. 

This model has a global Z 2 x Z 2 symmetry generated by 
the symmetry operations fJi even X > and IL odd X *- We 
consider this symmetry to be on-site, which is to say it 
acts on states as a unitary representation U(g) of the 
symmetry group G = Z 2 x Z 2 = {l,x,y,z} (with y = 
xz), such that U(g) acts as U(g) = [u(g)]® N , where N is 
the number of sites (we group qubits into two-qubit sites 
in order to ensure this condition is satisfied; see Fig. [I]). 
As we will see in Sec. |IIB| the cluster model lies in a 
nontrivial SPT phase with respect to this symmetry, so 
that the cluster state cannot be smoothly deformed into 
a product state without breaking the symmetry [26j . 

The 1-D cluster state can be represented as a pFCS 
[3J. For our purposes we will take the MPS tensor Ac to 
have the form 



Ac [++]=!, 
A C [-+] = Z, 



A c [+- 

M— 



] = x, 
■1 = xz 



-iY. 



(3) 



This is expressed with respect to a particular basis for 
a two-qubit site, where |±) = (|0) ± |l))/\/2. Here, and 
throughout this paper, we use the notation A[ip] to refer 
to the linear operator obtained from the MPS tensor A 
by interpreting 
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B. Symmetry-protected topological order in 
finitely-correlated states 

Here, we will review the results of [El [16] on the man- 
ifestation of SPT order in pFCS, and demonstrate that 
the 1-D cluster model indeed lies in a nontrivial SPT 
phase with respect to the Z 2 x Z 2 symmetry. 

Consider some ground state which is invariant under 
the on-site representation U(g) = [u(g)]® N of some sym- 
metry group G, and which can be represented as a pFCS, 
as m Eq. Q. The tensor A can be taken to obey a sym- 
metry condition [THl [55] 



A[u{g)^m=f3{g)W{g)^A{^}W(g), 



(5) 



where /3(g) is a one-dimensional linear representation of 
the symmetry group G, and W(g) is a projective unitary 
representation of the symmetry group G. This means 
that 



W(gi)W(g 2 ) = uj(gi,g 2 )W(g!g 2 ) 



(6) 



for some function lo, called the factor system of the pro- 
jective representation, which maps pairs of group ele- 
ments to phase factors. By blocking sites, we can ensure 
that (3(g) = 1 (however, for simplicity we will assume 
that (3(g) = 1 without blocking). Eq. ^ can then be 
represented graphically as 
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= - W(g) - A - W(g)t 



(7) 



u (g) 
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Observe that W(g) can be multiplied by a g-dependent 
phase factor without affecting Eq. ([7]); a set of factor 
systems related by such a transformation is referred to 
as a cohomology class. The arguments of [151 116j show 
that two such pFCS ground states correspond to the 
same cohomology class if and only if they are in the 
same symmetry-protected phase. Nontrivial cohomol- 
ogy classes [those not containing the trivial factor system 
u> (31,32) = 1] correspond to phases with nontrivial SPT 
order. 

As an example, consider the cluster model, and its 
Z2 x Z2 symmetry. The on-site representation u(g) of 
the symmetry is generated by 



u(x) 
u(z) 



X> 



X. 



(8) 
(9) 



and the MPS tensor is given by Eq. (|]). It can be shown 
that the symmetry condition Eq. Q is satisfied with the 
projective representation W — Vp, where Vp is the Pauli 
representation 



Vp(l)=I, Vp(x)=X, V P (z) 



, V v (y) = Y. 

( 10 ) 

This projective representation has nontrivial cohomol- 
ogy class, so that the cluster model lies in a non-trivial 
symmetry-protected phase. 



Symmetry-respecting perturbations to the 
cluster state 



Suppose we now consider a perturbation to the clus- 
ter Hamiltonian Eq. ([2]), such that the perturbed model 
still respects the Z2 x Z2 symmetry and admits a pFCS 
ground state. Unless the perturbation is large enough to 
induce a phase transition, the MPS tensor A correspond- 
ing to the perturbed pFCS ground state should still sat- 
isfy the symmetry condition Eq. (j7|, for some projective 
representation W(g) with the same factor system as the 
Pauli projective representation V [Eq. (10)]. 

The general form of the MPS satisfying these symme- 
try conditions was established in 13 . Here we will briefly 
review the relevant results from [T3] . We observe that the 
Pauli representation satisfies a property which we refer 
to as maximal non-commutativity: 



Definition 1. A projective representation W(g) of an 
abelian group G is called maximally non- commutative if 
the subgroup Z\y(G) = {g £ G : W(g) commutes with 
W(h) for all h £ G} (which we can think of as the "pro- 
jective centre" of G) is the trivial subgroup. 

Notice that the subgroup Z\y(G) is actually deter- 
mined by the factor system u>, since W(g)W(h) = 
Lu(g,h)W{gh) = ufehJu^sj-MW- Further- 
more, it is the same for all factor systems within a given 
cohomology class. Much of the discussion in this pa- 
per can be applied to any SPT phase characterized by 



a finite abelian symmetry group and a maximally non- 
commutative cohomology class. 

An important consequence of maximal non- 
commutativity of a factor system is |13| 

Lemma 1. For each maximally non- commutative fac- 
tor system u of a finite abelian group G, there exists a 
unique (up to unitary equivalence) irreducible projective 
representation with factor system ui. The dimension of 
this irreducible representation is \/\G\ . 

Specifically, the Pauli representation Vp of Z2 x Z2 is the 
unique irreducible projective representation correspond- 
ing to its factor system. In general, throughout this pa- 
per, we will use V(g) to denote the unique irreducible 
projective representation for the current factor system. 
A consequence of Lemma [T] is that, for a tensor satis- 
fying the symmetry condition Eq. Q, the bond space 
decomposes as a tensor product of a yj \ G | -dimensional 
protected subsystem in which W(g) acts irreducibly as 
V(g) and a junk subsystem in which W(g) acts trivially, 
i.e. 



W{g) = V{g)®l 



(11) 



Thus the tensor A appearing the MPS representation of 
ground states in the symmetry-protected phase satisfies 
the symmetry condition 



.4 



V(g) 



U Uv(,)tU 



.4 



(12) 



u (g) 



Here we use a thick line ( ) to represent the protected 

subsystem, and a dashed line ( — ) to represent the junk 
subsystem. 

Now suppose we perform a projective measurement 
on one site in a simultaneous eigenbasis {\i)} (which is 
|++), |H — ), I — h), I ) for the Z 2 x Z 2 cluster state sym- 
metry), and obtain the outcome Then the resulting 
state on the remaining sites is found by replacing the 
original MPS tensor at the measured site by 



.4 



(13) 



Now we make use of another consequence of maximal 
non-commutativity [13] : 

Lemma 2. Let u(g) be a linear on-site representation 
of a finite abelian symmetry group G, and let to be a 
maximally non- commutative factor system of G. Then 
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for each basis element \i) in a simultaneous eigenbasis 
{\i}} ofu(g), there exists a group element gi such that 



Xi(9)V(g) = V( 9l )V(g)V( gi 



VgeG 



(14) 



for any projective representation V(g) with factor system 
u>, where Xi(d) * 5 ^ e scalar representation of G such that 

u{g)\i) = Xi(g)\i)- 

For the example of the cluster state symmetry, we have 

<7++ = 1, <7_| = x, g |- = z, g = y, as can readily be 

verified directly. 

As was shown in |13j . Lemma[2]in conjunction with the 
symmetry condition Eq. ( |12[ ) implies the decomposition 
A[i] = V{gi) ® A[i], represented graphically as 
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VQh) 



A 



(15) 



for some tensor A. Another way of writing this result is 
that 



.4 



(16) 



where we have defined the tensor 

v y— ( A v{ 9i ) V\ 



E 



(17) 



/ 



Note that, from a quantum circuit perspective, this ten- 
sor can also be interpreted as a unitary controlled opera- 
tion J2i \i)(i\®V(9i) coupling a site to an ancilla particle; 
hence the choice of notation. Conversely, any MPS ten- 
sor of the form Eq. (16) for some tensor A will satisfy 
the symmetry condition Eq. (12). Following [25], we re- 



fer to the tensor A as as the degeneracy tensor; and to 
the tensor of Eq. 



the symmetry, as t 

From Eq 
picture 



17), which is determined entirely by 



(15) 



le structural tensor. 
we see that, in the correlation space 



measuring in a simultaneous eigenbasis {\i}} 
leads to an identity evolution (up to outcome-dependent 
byproduct operators) on the protected subsystem of cor- 
relation space; this evolution is determined by the sym- 
metry (hence the same throughout the SPT phase), and 
decoupled from the junk subsystem. Therefore, we say 



that the identity gate operates perfectly throughout the 
SPT phase. However, the result of measurement in 
any other basis is not fixed by the symmetry, and in 
general leads to the protected subsystem being coupled 
to the junk subsystem, so that the operation of other 
measurement-based gates will not be a robust property 
of the symmetry-protected phase. 



D. The dual picture for MBQC on a ID resource 
state 

In order to deal with the randomness of measurement 
outcomes, the measurement protocol for MBQC with the 
cluster state needs to be adaptive: the outcome of the 
measurement on one site will affect the measurement ba- 
sis on other sites arbitrarily far away. In analysing the 
effect of this protocol when acting on a perturbed resource 
state, we would like to make an argument based on the 
locality of the perturbed Hamiltonian, but the non-local 
adaptivity of the measurement protocol poses a difficulty. 
Therefore, in this section, we develop an alternate char- 
acterization of the effect of the cluster state adaptive 
measurement protocol acting on a ground state in the 
symmetry-protected phase. We will show that this pro- 
tocol is equivalent to a dual process acting on a related 
state, which we call the dual state. We will show that this 
dual process simply consists of a sequence of unitary in- 
teractions between selected sites (those corresponding to 
the locations of non-trivial gates) and an ancilla particle, 
with no adaptivity. 

In our discussion of the dual process, we will represent 
a pFCS ground state on an infinite chain as a formal 
tensor network 



I*) = 
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This is not, of course, the mathematically rigorous way to 
describe pFCS, but we find it useful for facilitating under- 
standing. In Appendix [X] we will discuss how to formu- 
late similar arguments in the rigorous pFCS framework. 
Later on (in Section [nT| , we will also be interested in fi- 
nite chains; in that case, the arguments of this subsection 
can be applied more directly, given appropriate bound- 
ary conditions [specifically, the boundary conditions at 
the right edge should be as depicted in Eq. (51)]. 

The usefulness of the exact 1-D cluster state (with 
MPS tensor Aq) as a quantum computational wire re- 
sults from the fact that, for each gate U in a certain set, 
there exists a basis {|a)} such that 



— A c — = — U —B a — 



(19) 



where B a 
operator. 



is the outcome-dependent unitary byproduct 
When we measure one site projectively and 
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obtain the outcome |a), the original MPS tensor A is 
replaced at the measured site by Eq. ( 19 ) in the tensor- 



network description of the resultant state. 

In the case of the exact cluster state, the effect of the 
byproduct operator can be accounted for by adjusting 
the measurement basis for future measurements. This 
fact turns out to be closely related to the nontrivial SPT 
order, as we now demonstrate. Our discussion relies on 
the observation that, in the cluster state, the byproduct 
operators are Pauli operators. That is to say, it is always 
the case that B a is a scalar multiple of V(g a ) for some 
9a € Z2 x Z2. Hence, we can make use of the symmetry 
condition [which can be derived from Eq. ( 12 )] 



— B a — Aq — — Aq — B a — 



(20) 



b a 



where b a = u(g a ). Applying this condition repeatedly 
shows that the byproduct operator can be displaced ar- 
bitrarily far to the right. In our formal tensor-network 
picture for an infinite chain, we consider that this pro- 
cess is continued indefinitely, so that the byproduct op- 
erator "disappears out to infinity" , and is replaced with 
b a acting on all sites to the right of the one on which the 
measurement took place, i.e. 

— s^Wa^a^-Ia^Wa^— • • • 



€E G. Thus, the resource state is of the form 



.4 



I*) 



.4 
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.4 



.4 



(22) 

with the MPS tensor A of the form Eq. (|16|). We now 



repeat the above argument, in reverse. We make use of 
the symmetry condition Eq. ( 12 ) in the form 



.4 



.4 



(23) 



bi 



from which we obtain 



.4 







.4 



.4 



.4 



.4 



6* &t 6* 



Ac— An— Ac— A 



(21) 



b a 



b a 



b a 



Hence, whenever we obtain the "wrong" outcome for a 
measurement (i.e. the corresponding byproduct opera- 
tor B a is not the identity), we can recover the "correct" 
resultant state by applying the correction to all the 
remaining sites on the right (equivalently, we can simply 
adjust the measurement basis for measurements on those 
sites). 

Let us now examine what happens when we perform 
the same adaptive measurement protocol on a resource 
state that is not the exact cluster state. In general, we 
can consider any pFCS ground state contained within 
any SPT phase characterized by finite abelian symmetry 
group G and a maximally non-commutative cohomology 
class, but we assume that the adaptive correction takes 
the same form as for the cluster state, i.e. application of 
u(g a ) to the sites on the right for some group elements 



■Bt 



A 



A 



A 



A 



A 



(24) 



Therefore, we have shown that the process we actu- 
ally perform, i.e. applying the measurement-dependent 
correction to the sites on the right of the one mea- 
sured, is equivalent to a different process, in which the 
measurement-dependent correction is applied in the in- 
ternal bond space of the MPS, as depicted in the right- 
hand side of Eq. ( 24 ) . In a physical system, of course, we 



do not have direct access to the internal bonds of a ten- 
sor network state, so we could never perform the latter 
process directly; nevertheless, the two are equivalent. 

Following the measurement and the adaptive correc- 
tion, which we think of as being performed internally, as 
in the right-hand side of Eq. (24), the outcome of the 



measurement can be "forgotten" , i.e. we describe the re- 
sultant state of the system as the mixture of the right- 
hand side of Eq. ( 24 ) for all possible measurement out- 
comes. Without affecting the reduced state on the re- 
maining unmeasured sites, for notational convenience we 
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replace this mixture with a coherent superposition, i.e. 





-6ft- 
















A 






A 




A 




A 




A 



(25) 



where we have defined the tensor 



— v 

P 



E 



[a I 



(26) 



/ 



which we can also interpret as a unitary coupling 
Y, a \a){a\ ®Bj = X) a l a )H ® ^(ff") f coupling a site to 
an ancilla particle. We now use the k index throughout 
the paper to distinguish the "G" 's resulting from different 
measurement operations. (The label k refers to the site 
at which the measurement is being performed; we include 
this label to reflect the dependence on the measurement 
basis {la)} and byproduct operators B a , which will in 
general be different for each site at which a measurement 
is performed.) 

Now, using the expression Eq. ( 16 1 for the MPS tensor 
A, we can write 







A 






I 







(27) 




(b) 



A - A -- A -- A -- A 



— 



FIG. 2. (a) After the adaptive measurement sequence, we can 
treat the resultant state as having the form shown. This is 
equivalent (b) to building |$) from the dual state |^) [the 
shaded box; see Eq. ( 29 1] by unitary couplings to an ancilla 
particle. 



the state built from the degeneracy tensor A, with the 
structural tensor discarded: 



I*) 

















A 




A 




A 




A - 



(29) 



where we have defined the tensor Gk (which can also be 
interpreted as a unitary coupling between a site and the 
ancilla particle) according to 




-®-@- 

* A- 



(28) 



We are now in a position to define our dual process. 
Suppose we perform a sequence of such adaptive mea- 
surements at successive sites, which at each site is de- 



scribed by the insertion of the tensor Eq. (26), as in Eq 
(|25l 



There will be a different coupling Gk associated 
with each site k. As shown in Figure [2j we find that the 
original adaptive measurement process, applied to the re- 
source state |\f r ), is equivalent to a dual process applied 
to the dual state |^). The dual state |\1/) is defined to be 



The dual process comprises a series of consecutive unitary 
interactions Gk between individual sites k and an ancilla 
particle. 

There are several reasons why this "dual picture" is 
a useful way to understand the operation of MBQC in 
one-dimensional ground states. First, the dual process 
lacks the long-range measurement adaptivity which is a 
characteristic of the original adaptive measurement pro- 
tocol. Second, the perfect operation of the identity gate 
is automatically built in, because, for sites k at which the 
adaptive measurement process at the given site is the one 
that corresponds in the exact cluster state to the iden- 
tity gate, [i.e. the measurement basis is the simultaneous 
eigenbasis {\i)} of the symmetry, and the byproduct op- 
erators are Bi — V(gi), where the gi are the group ele- 
ments appearing in Eq. ( 15 )], the corresponding coupling 
is trivial, Gk = !• 

The final motivation for the dual picture is that the 
dual state on which it is based has some physical signifi- 
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cance in its own right, and retains some key properties of 
the the original resource state. For example, if the orig- 
inal MPS tensor A generates a pFCS, then so does A, 
and the respective correlation lengths obey the inequal- 
ity £ < £ (see Appendix |a| . Additionally, in Appendix 
|B]we will show how our dual state can be obtained from 
the original ground state through a generalization of the 
unitary that was introduced by Kennedy and Tasaki [30] 
to transform the SPT Haldane phase [3TJ [32] into a lo- 
cal symmetry-breaking phase; this unitary has recently 
been described as a "topological disentangler" [33], and 
in some sense we can think of the dual state as being a 
topologically disentangled version of the original resource 
state. 



The 



process for initialization and readout in the 
1-D cluster state 



Above we only discussed measurement sequences cor- 
responding to unitary gates in correlation space. A com- 
plete scheme for using a 1-D resource state as a quantum 
computational wire also includes measurement sequences 
corresponding to initialization (i.e. discarding the cur- 
rent state in correlation space and replacing it with a 
fixed state), and readout (i.e. making the state in corre- 
lation space available as the physical state of one qubit). 
We now describe briefly how the measurement protocols 
used on the 1-D cluster state for these purposes can be 
accommodated in our framework. 

Initialization. — The initialization procedure in- 
volves measuring a site in the computational basis 
{|00), |01), 1 10) , |11)}. In this basis, the MPS tensor A c 
for the exact cluster state takes the form 



(30) 
(31) 
(32) 
(33) 



Ac [00] = 
A c [01] = |1)(0|, 
A c [10] = |0)(1|, 
A C [11] = -|1)<1| 



The randomness of measurement outcomes can therefore 
be accounted for by applying the appropriate outcome- 
dependent correction operator in correlation space fol- 
lowing the measurement: £>q = b\ q = I, = B\ x = 
X. Since the correction operators are Pauli operators, 
the above discussion applies without change. 

Readout. — The standard readout procedure for the 
cluster state involves measuring the second qubit of a 
two-qubit site in the computational basis, then applying 
an outcome-dependent correction operator to the first 
qubit, which acts as the output. Provided that we are 
only interested in the final state of the output qubit, this 
procedure is equivalent to a coherent correction operator 
coupling the two qubits in the site (specifically, it is a 
controlled-Z gate C z = |0> <0j ®I+ |1)(1| <g) Z). Carrying 
through a similar argument to that given above for uni- 




the interaction Qk in the dual process between the site k 
in question and the ancilla particle given by 



(34) 



E. MPS of minimal bond dimension and the dual 
picture 



As an example of the general formalism introduced 
in Sec. II D here we will examine the form of the cou- 




plings Qk appearing in the dual process [Eq. ( 28 )] , in 



the particular case that the res ource state is an MPS 
with bond dimension D = ->J\G\, where G is the symme- 
try group characterizing the symmetry-protected phase. 
Given that the dimension of the protected subsystem is 
\f\G\ (by Lemmall|, this is the smallest possible value of 
D, and corresponds to the absence of a junk subsystem 
(or, more precisely, a junk subsystem of dimension 1). In 
particular, the ID cluster state is of this type. In general, 
the MPS tensor A for such an MPS must be of the form 



A\i] = A[i]V{ 9i 



(35) 



where the A[i] here are scalars. It follows that the dual 
of such a state is a product state, with each site in the 
state \4>) = Yli A-WN)- (We choose the normalization for 
the MPS tensor A so that {<j>\<j)) = 1.) Therefore, the 
effect of the dual process acting on the dual state results 
from a series of independent interactions of the form 



\<t>) 



(36) 



We recall that, in the correlation space picture of quan- 
tum computational wires, a quantum state can serve as 
a resource for executing a unitary gate U if there exists 
some basis {\a}} such that 



A[a] = p a B a U, 



(37) 



for some set of unitary byproduct operators B a and 
scalars j3 a . We will now show how this property man- 
ifests itself in the dual picture, for the class of states 
considered here. We make use of the representation for 
the MPS tensor A as 



tary gates, Sec. II D , we obtain the same result, but with 



(38) 



It follows that, at a site k measured in the basis {|a)}, 
with the byproduct operators B a , we have 



(39) 





U 



(40) 
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where \<p) = J2 a Pa\&). (It can be shown that our choice 
of normalization ensures that (tp\tp) — 1.) Here the first 
equality follows from Eq. (38) and the definition of Q; 

Thus, we 



and the second inequality follows by Eq. ( |37[ ) 
have shown that in the dual picture the gate U simply 
acts on the ancilla particle. 

Next we will do a similar analysis for the initialization 
and readout procedures specific to the 1-D cluster state. 

Initialization. — Recall the discussion of initialization 
in Sec. Ill Dl We make use of the form of the MPS ten- 
sor Ac in the computational basis, Eqs. (30-33), mul- 



tiplied by the appropriate normalization factor as dis- 
cussed above. Thus, for the site k at which initialization 
takes place, we find that 



(41) 





(42) 



where T = ^ (|00) (0| + |01) (0| + |10) (1| - |11) (1|) 
(thanks to our choice of normalization, we find that T is 
an isometry, i.e. T^T = I). Therefore, applying the mea- 
surement sequence for initialization leads to the ancilla 
system getting initialized in the state |0), as we would 
expect. 

Readout. — From the definition of the operator Qk in 
the case of sites k at which readout takes place [Eq. (34)], 



V^fc 















J 
















1*) 














J 















»5 




(a) 



FIG. 3. (a) As long as the locations of non-trivial gates (sites 
k and I in this diagram) are separated by a distance much 
greater than the correlation length f, the reduced state on 
those locations will be approximately a product state, and the 
dual process then reduces (b) to a sequence of independent 
interactions. 



stituent qubits, each denoted by a thick line) 




■Ac 



:c z 



(43) 



(44) 



we find that (here we separate a site into its two con- 



where | J) = (l/S/2) (|00) + |11» is the canonical 
maximally-entangled state. Thus the state of the ancilla 
qubit indeed gets transferred onto the output qubit. 



F. MBQC on a perturbed resource state simulates 
a noisy quantum circuit 

In the previous subsection, we saw how measurements 
on an MPS of minimal bond dimension correspond to 
quantum gates. Now we will consider what happens when 
we perform the same measurement sequences on a per- 
turbed resource state, assuming that the perturbed state 
remains within the same SPT phase. We will find that 
measurements on such a perturbed cluster state simulate 
the same quantum circuit, but with noisy gates. The 
noise is described by application of a completely positive, 
trace preserving (CPTP) noise superoperator following 
each gate. 

In Section [II E[ we were able to treat each gate inde- 
pendently in the case of the unperturbed cluster state 
because the dual state was a product state, |\?) = \</>)® N . 
This will no longer be true once we introduce pertur- 
bations, but we still want to treat gates independently 
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Towards this end, we recall that for a site for which 
the corresponding sequence is that for the identity gate, 
the associated coupling Q k in the dual picture between 
that site and the ancilla particle is trivial. Therefore, 
such a site can be traced out from the beginning with- 
out affecting the final state of the output qubit. That 
is to say, we only need to consider the reduced state 
p = Tr tr iviai sites |*) (^1 on the remaining sites, which 
are those corresponding to non-identity gates (we refer 
to these as the non-trivial sites). We are free to choose 
our measurement protocol to ensure that the distance 
between any two non-trivial sites is much greater than 
the correlation length £. In that case, we expect that 
the reduced density operator p should be approximately 
a product state between the non-trivial sites 



Pk — Pprod- 



(45) 



k non— trivial 



Indeed, for pFCS, it is straightforward to show that there 
exists a choice of p pro d such that 



\P - Pprodill < mf(R), 



(46) 



where ||.||i denotes the trace norm, m is the number of 
non-trivial sites, R is the minimum distance between any 
two non-trivial sites, and f(R) is a function related to the 
transfer channel of the pFCS, behaving asymptotically 
as f(R) = 0[exp(—R/£)] for large R, where £ is the 
correlation length associated with the pFCS. 

We first consider the case where p = p pro d exactly. 
Because p is then a product state, we find, as in the 
previous subsection, that the dual process acting on the 
dual state is again effectively a sequence of independent 
interactions, this time of the form 



Pk 



(47) 



(see Fig. [3|. Let us suppose that Qk results from the 
measurement sequence corresponding to a unitary gate 
Uk- Then, after tracing out the physical site, Eq. (47) 



corresponds to an evolution on the ancilla qubit described 
by the CPTP map 



Ak{<j) = Tr physical site [Gk{o- <g) Pk)G\ 



(48) 



As we saw in Sec. |II E[ in the absence of perturbations 
to the cluster state, pk — I^K^I an d <Afc = Uk, where 
Uk(<j) = UkcrUl. In general we can write Ak = £k °Uk, 
where £ k is a noise superoperator for which it is straight- 
forward to show that 



\\£ k 



o 



< 



\Pk 



(49) 



where ||.||<> is the diamond norm on superoperators |34j . 
The cases when Qk corresponds to initialization or read- 
out are analogous. Therefore we have shown (in the case 



P = Pprod) that the measurement protocol on the per- 
turbed cluster state reproduces the desired quantum cir- 
cuit, except that each gate (as well as the initialization 
and readout steps) is accompanied by some associated 
noise. Furthermore, if the perturbation is sufficiently 
small then we expect that each pk will be close to \(f))(<f)\, 
so that the noise will be weak, in the sense that £k is 
close to the identity superoperator in the diamond norm. 

In the general case, in which p and p pro d are not equal, 
but are e-close in the trace distance, we just need to ob- 
serve that the reduced state of the output qubit following 
the dual process can be obtained from p by application of 
some CPTP superoperator, which we call B. From the 
contractivity property of the trace distance, it follows 
that \\B(p) - B(p P rod)||i < \\p- /Qprodlli < £■ Therefore, 
the effective noisy quantum circuit description correctly 
describes the final state of the output qubit up to an ac- 



curacy e. Note that, because the bound Eq. (46) depends 
on the number of non-trivial gates m, it will be necessary 
to have the separation R scale with m in order to obtain 
a fixed accuracy e, but only logarithmically; indeed, the 
minimum separation required to achieve an accuracy e 
scales like -R m in/£ = C[l°g(m/e)]. 



G. Summary of Section [IT] 

In Section [TTJ, we have presented, within the context 
of pure finitely-correlated states, the main ideas leading 
to our effective noise model construction. Our discus- 
sion has hinged around the "dual state" which we asso- 
ciated with each ground state carrying the appropriate 
SPT order. Loosely speaking, we can think of the en- 
tanglement in SPT-ordered ground states as comprising 
"topological" and "non-topological" components inter- 
twined. The topological component is fixed throughout 
the phase and is responsible for the distinctive charac- 
teristics of the SPT phase, such as the degeneracy in the 
entanglement spectrum |35j , the diverging localizable en- 
tanglement length |3H [37] , and the perfect operation of 
the identity gate. One can think of the dual state as being 
obtained from the original ground state by a topological 
disentangler, "separating out" the topological component 
of the entanglement and leaving only the non-topological 
component |33j . 

In this paper, the importance of the dual state is due 



to the following fact, which we established in Sec. II D 



the cluster state adaptive measurement protocol, when 
applied to an SPT-ordered ground state, couples in a 
natural way to the topological component of the entan- 
glement, and the effect is thus equivalent to a "dual pro- 
cess" (with a simpler structure) acting on the dual state. 
This result gives rise to an effective quantum circuit de- 
scription describing the outcome of the measurement pro- 
tocol applied to any SPT-ordered ground state [e.g. see 
Fig.^a)]. The action of non-trivial gates is determined 
by an interaction with a single site in the dual state, and 
perturbations to the dual state give rise to noisy gates. 
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As long as the locations of non-trivial gates are suffi- 
ciently separated, the reduced state on the sites relevant 
for the gate operation will be a product state, and this 
corresponds to independent noise acting on each gate in 
the effective circuit description. 



III. THE EFFECTIVE NOISE MODEL 
CONSTRUCTION: GENERAL GROUND STATES 

In this section, we will extend our characterization of 
the effective noise model to any ground state within the 
symmetry-protected phase, without reference to finitely 
correlated states. Instead of starting from scratch, 
we will build on the results of Section [Hj as follows. 
We formulate a condition which we believe (on physi- 
cal grounds) to be satisfied for any system within the 
symmetry-protected phase. We will show that this con- 
dition leads to a construction for the dual state of any 
ground state in the SPT phase, independently of the 
pFCS formalism. Furthermore, given an exact MPS rep- 
resentation for the dual state (which always exists, albeit 
possibly with a bond dimension exponentially large in 
the system size), we show that one can construct a cor- 
responding MPS representation for the original ground 
state, such that the arguments of Section [TT] can be ap- 
plied without significant change. In order to establish the 
approximate factorization condition Eq. ( 45 ) in the case 



of general ground states, we will show that the dual state 
is (like the original resource state) the gapped ground 
state of a local Hamiltonian, which can be constructed 
in a straightforward way from the original Hamiltonian. 
This will allow us to establish the approximate factoriza- 



tion condition Eq. ( 45 1 without assuming that the dual 



state has a pFCS structure. 



A. Symmetry-protected topological order and 
boundary conditions 

Because we are considering general ground states, we 
can no longer make direct use of the characterization of 
SPT order in finitely-correlated states of [THl HE|- In- 
stead, we adopt the perspective in which SPT order is 
related to the fractionalized edge modes associated with 
open boundary conditions |32| . Our discussion will, out 
of necessity, be physically motivated rather than math- 
ematically rigorous, but will suggest the formulation of 
the precise assumptions under which the rigorous results 
of this paper can be proven. 

Consider a 1-D chain with open boundary conditions, 
with symmetry-respecting interactions such that, in the 
bulk, there is no symmetry-breaking and a finite en- 
ergy gap for excitations. It is still possible that the 
energy gap for edge excitations may be much smaller 
than the bulk gap (or even zero). We denote by V the 
subspace comprising the low-lying edge states. If the 
chain is sufficiently long, we expect that the gap in the 



bulk should ensure that the edges are non-interacting, so 
that V decomposes a tensor product of degrees of free- 
dom associated with the left and right edges respectively, 

V = Vi ®V r , and the restriction H-p of the Hamiltonian 
H onto this subspace is a sum of terms acting on each 
edge: H v = hi ® I + I ® h r . 

For our purposes, we are interested in how the edge 
states transform under the symmetry. Because V is 
spanned by a set of energy eigenspaces, and the Hamilto- 
nian commutes with the representation of the symmetry 
U(g), it follows that V must be an invariant subspace 
for U(g). We write the operation of U(g) restricted to 

V as U-p(g). We expect that, for sufficiently long chains, 
the symmetry should act independently on the respective 
edge states, i.e. 



U v (g) = V l (g)®V r (g). 



(50) 



By assumption, U{g), and hence U-p(g), is a linear repre- 
sentation of the symmetry group G. It follows that Vi(g) 
and V r (g) must be projective representations of G, and 
if Vi(g) has factor system uj, then V r (g) must have factor 
systemu; -1 . We are free to transform V r (<?) — » f3(g)V r (g), 
Vi(g) ~^ Z^ 1 (5)^(5) f° r an Y g-dependent phase factors 
(3(g) without affecting Eq. (50), but the cohomology class 
[uj] is uniquely determined. Furthermore, it is intuitively 
clear that any continuous symmetry-respecting variation 
in the Hamiltonian cannot change the cohomology class 
[uj], except at a phase transition (where the gap closes 
in the bulk, and the left- and right- edge modes need no 
longer be non-interacting). Therefore, we have an alter- 
native characterization of the SPT phase corresponding 
to a cohomology class [uj] : it comprises the systems where 
a left edge is associated with emergent edge states trans- 
forming projectively under the symmetry with cohomol- 
ogy class [uj]. This is a generalization of the well-known 
observation that systems in the S l O(3)-invariant Haldane 
phase have emergent spin- 1/2 degrees of freedom at the 
edges [5E1GH]. 

In non-trivial SPT phases, the edge interactions hi 
and h r (and therefore the overall Hamiltonian H) will 
always have degenerate ground states, due to the fact 
that non-trivial projective representations cannot be one- 
dimensional. On the other hand, we expect, at least in 
the case that the symmetry group G is abelian, that a 
non-degenerate ground state can be recovered by intro- 
ducing terminating particles at the left and right edges, 
transforming projectively under the symmetry with fac- 
tor systems uj' 1 and u respectively (see Figure [3]). This 
is because, loosely speaking, these terminating particles 
can couple to the edge modes, with the composite system 
at each edge transforming under a linear representation 
(and therefore, in the case of an abelian symmetry group, 
generically having a non-degenerate ground state). For 
example, the ground state of a spin chain in the Haldane 
phase can be made non-degenerate through coupling to 
spin-1/2 particles at the edges. Conversely, if the termi- 
nating particles do not transform with the cohomology 
classes [uj' 1 ] and [uj] respectively, then the degeneracy 
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(a) 



m z i 

(a) Open boundary conditions 



LO id UJ CJ 

(b) With terminating particles 

FIG. 4. (a) The low-lying energy subspace V of a 1-D chain 
with open boundary conditions, in an SPT phase character- 
ized by the cohomology class [u] , decomposes as a tensor prod- 
uct of "emergent edge modes" (red) associated with each end, 
transforming projectively under the symmetry, (b) The de- 
generacy can be removed by coupling terminating particles 
(green) at each end, leading to an effective coupling to the 
edge modes. 



cannot be removed completely because there is still a 
non-trivial projective symmetry transformation at each 
edge. 

Thus, we have arrived at yet another characterization 
of SPT order, which we state as a conjecture in the ab- 
sence of a rigorous proof: 

Conjecture 1. A 1-D chain respecting an on-site repre- 
sentation of an abelian symmetry group G is in the SPT 
phase characterized by cohomology class [to] if and only 
if the following condition is satisfied: 

Condition 1. The finite-chain ground state can be 
made non-degenerate and gapped by the inclusion of 
symmetry-respecting interactions coupling the left and 
right edges of the chain to terminating particles trans- 
forming projectively under the symmetry, with factor 
systems ui and w _1 respectively. 

In any case, in the remainder of this section, we will 
consider systems satisfying Condition [I] Specifically, all 
the results will apply to finite chains with the appropri- 
ate edge couplings imposed to ensure a non-degenerate 
gapped ground state. This will prove convenient for our 
analysis, but the properties of the system in the bulk 
should not, of course, depend on the boundary condi- 
tions. 

Note also that, in the case of a system with the interac- 
tions governed by the parent Hamiltonian of a pFCS [gen- 
erated by an MPS tensor satisfying the symmetry con- 
dition Eq. ([7| corresponding to the symmetry-protected 
phase] , Condition [l] can easily be established directly. 
Furthermore, the stability theorem of [40] ensures that 
Condition[l]rcmains true for sufficiently small symmetry- 
respecting perturbations of such models, regardless of the 
validity of Conjecture [I] 



YD 



(b) 



0-0- -{a} -{a} -{a} -\a}[r_ 



FIG. 5. The "topological disentangler" T> applied to the orig- 
inal ground state (a) turns it into the dual state (b), leaving 
the terminating particles maximally entangled. 



B. The general construction for the dual state; 
exact MPS representation of SPT-ordered ground 
states 



Recall that in Section II D we defined the dual state 



in the context of pFCS. Here, we will give an analogous 
construction for the dual state corresponding to a general 
ground state within a symmetry-protected phase, pro- 
vided that the phase is characterized by a finite abelian 
group G and a maximally non-commutative cohomology 
class [uj] . The construction applies to a finite chain, with 
the appropriate boundary conditions as discussed in Sec- 
tion III A This construction will then allow us to express 



the original ground state as an MPS, with the MPS ten- 
sors satisfying an appropriate symmetry condition. 

We consider a finite chain coupled to terminating par- 
ticles, such that the overall system is invariant under the 
symmetry U(g) = V*(g) <g> [u(g)]® N <g> V(g). Here we 
have taken the right terminating particle to transform 
under V(g), the unique irreducible projective representa- 
tion with factor system lo; and the left terminating par- 
ticle under V* (g) [V*(g) is the operator obtained from 
V(g) by complex conjugation of the matrix elements in 
some basis; observe that V*(g) is a projective represen- 
tation of G with factor system cj~ ]. 

The natural analogues in the current setting (finite 
chains, with the specific choice of boundary conditions) 
of the pFCS ground states which we considered in Section 
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are states of the form 



I*) 







.4 



T 



A 



A 



A 



A 



R 



(51) 



for some end vectors \L) and (R\, and where the MPS 
tensor A satisfies the symmetry condition Eq. (l2| [which 
can be shown to ensure the invariance of the state under 
[/(<?)]. Given the decomposition Eq. (161 for the MPS 



tensor A, it follows that the dual state can be obtained 
from the original ground state by a sequence of unitary 
interactions between individual sites and the terminat- 
ing particle on the right (see Fig. [5]); we can think of the 
overall unitary transformation T> as a "topological disen- 
tangler". Specifically, we have T>\^>) = |\&) ® |7), where 

1 7) = X)fe=i 1^)1^) i s the- canonical maximally-entangled 
state between the terminating particles. 

We will now show that, for a general gapped symmetry- 
respecting ground state |\&) [not necessarily in the MPS 
form Eq. (51 1], it remains the case that X>|\P) = |$) <& |7) 
for some state [*) on the non-terminating sites; this will 
serve as the definition of the dual state |^/) for general 
ground states. 

We observe that the original ground state | \&) must be 
invariant under the global symmetry operation U(g), i.e. 



U(g)\*)=a(g)\V). 



(52) 



(Without loss of generality, we can set a (g) = 1 by ab- 
sorbing it into into the definition of the symmetrjf^]) This 
implies that the state X>|\E r ) is invariant under T>U{g)T>' . 
Let us examine what this 'dual' symmetry looks like. We 
observe that 

V 1 [u{g) ® V( 9 )]V\ = X)xi(ff)NX*l ® V(^V(g)V(g t ) 

(53) 

= J2\i)(i\®V(g) (54) 

i 

= I®V(g), (55) 

where 

(vt 



V l = 



(56) 



is the interaction from which T> is built; to get to Eq. 
(54), we made use of Eq. (14). From this, one can show 



that 



VU{g)V^ = V*(g) ® I® N ® V(g). 



(57) 



1 Specifically, we replace the action of the symmetry V(g) on the 
right terminating particle with its equivalent under a unitary 
transformation (by Lemma 1 of Ref. 13 ), a" 1 (g)V(g). 



It it straighforward to show (using the irreducibility of V) 
that invariance of a state under the right-hand side of Eq. 
(57) implies that it must be of the form V\^f) = |$) ® |7) 



for some state \ty), as required. 

It is now straightforward to construct an appropriate 
MPS representation for a general ground state. Indeed, 
let us consider an MPS representation for the dual state 
1$) of the form 



I*) 



A, 



A, 



Aa 



(58) 



We choose this representation to be exact; this may re- 
quire the bond dimension to be very large (scaling expo- 
nentially in the system size), but that is of no importance 
to us. Then we have 



2?|*) = |7)<g>|#) 



M 




A 2 




A 3 




Ai 




A 5 



(59) 



(60) 



Now we can apply the inverse transformation to ob- 
tain 



I*) 



<^ — &i — &i — & — ® 



A 2 



A 3 



(61) 

This is a representation of |W) as an MPS, with each of 
the shaded regions corresponding to an MPS tensor A 
of the form Eq. (16), and hence satisfying the symme- 



try condition corresponding to the symmetry-protected 
phase. In addition, we should take note of the boundary 
conditions at the right edge. These boundary conditions 
ensure that the arguments of Section |II D| apply without 
any need to invoke an infinite-system limit. 



The dual state as the ground state of a local 
Hamiltonian 



In the previous subsection, we have constructed the 
dual state for any ground state in the symmetry- 
protected phase. The original ground state is, by as- 
sumption, the gapped ground state of a local Hamilto- 
nian. In this subsection we will show that this is also 
true of the dual state. That is, starting from the original 
Hamiltonian 77, we construct another local Hamiltonian 
77 for which the dual state is the gapped ground state. 
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We start by proving a useful fact about the unitary 
transformation T> introduced in the previous section: al- 
though it is in general non-local, it maps symmetry- 
respecting local observables (i.e. those supported on a 
small set of sites of finite size) to local observables. In- 
deed, let us consider some local observable h; we will 
show that T>hT>' is also local. For concreteness, we sup- 
pose that h acts on two adjacent sites somewhere in the 
bulk. Now, observe that 




(62) 



where 




(63) 



By means of Eqs. (63) and (55), it can be verified that 



if h commutes with the symmetry, i.e. 







h 




h 




u (g) 


, (64) 














u (g) 





then 



v{g) 



V(g) 



(65) 



Since V(g) is an irreducible projective representation, Eq. 
(65) implies (by Schur's Lemma) that h acts trivially on 
the terminating particle, i.e. 



(66) 



Now, using Eqs. (62) and (66), we find that VhTfi = h, 



where h acts on the same two sites as h (see Fig.|6j). Thus, 
although the duality transformation D is non-local, we 



{yt) — (yt) — (yt) — (yt) — (yt) — (yt} 
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FIG. 6. From Eqs. (62 1 and (661, we get the pictured equal- 
ity. This shows that T>h — hT>, or equivalently DKD^ — h. 



have shown that it maps local symmetry-respecting op- 
erators to local operators, as promised. The exception 
is operators h at the left edge, which act non-trivially 
on the left terminating particle; in that case, the above 
argument breaks down, but we can observe directly from 
the structure of V that VKD^ is supported on the union 
of the support of h and the right terminating particle. 
For operators h acting non-trivially on the right termi- 
nating particle, the argument must be adjusted, but the 
conclusion that T>KD^ is supported on the support of h 
still holds. 

We are now in a position to construct the Hamilto- 
nian for which the dual state j*) is the gapped ground 
state. We observe that VHV' has |7) (g) as its gapped 
ground state; however, it includes terms acting non- 
trivially on the terminating particles. We define a Hamil- 
tonian acting only on the intermediate sites according to 
H = (I\VHV^\I) = T(H); by the locality result proven 
above, each local interaction term in H corresponds to 
a local term in H supported on the same set of sites. It 
can be shown that |*) is the unique ground state of H, 
and that the gap is at least as large as that of T)HT>' , or 
equivalently H. 
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D. The factorization condition for general ground 
states 

Recall that the other condition that needed to be sat- 
isfied in order to apply the arguments of Section [IT] for 
general ground states was that the factorization condi- 
tion for the reduced density operator p on the non-trivial 
sites in the dual state, 



) Pk = Pprod 



(67) 



should be satisfied when the non-trivial sites are suffi- 
ciently separated from each other. Recall that, for the 
case of pFCS, one can prove the bound 



||p - PprodHi < rnf(R), 



(68) 



with m the number of non-trivial sites, and f{R) a func- 
tion scaling asymptotically as f(R) — 0[exp(— R/£)], 
where £ is the correlation length in the dual state. We 
conjecture that Eq. (681 should be a general property of 



all gapped ground states of local Hamiltonians. However, 
we have only been able to rigorously prove the weaker 
bound 



||p-Pprod||i <md 2m f{R), 



(69) 



where f(R) is as before, and d is the dimension of the 
Hilbert space at each site; see Appendix [C] for the proof. 
Note that if we assume only the weaker bound Eq. 



(691, then the separation between non-trivial sites will 



need to scale more rapidly with the number of gates m; 
we find that the minimum separation i? m i n required for 
an accuracy e scales like 



R mi Jti = 0{m)+0[\og{l/e)]. 



(70) 



This still implies that the number of measurements need 
scale only as a polynomial in the number of non-trivial 
gates. 



E. Nonzero temperature 

The formulation of the dual state as the ground state of 
a dual Hamiltonian extends naturally to nonzero temper- 
ature: under the topological disentangler T>, the thermal 
state of the original Hamiltonian H maps to the thermal 
state of a dual Hamiltonian H. Furthermore, it can be 
shown that an appropriate adaptive measurement proto- 
col acting on the thermal state of the original Hamilto- 
nian is equivalent to a non-adaptive dual process (of the 
same form as in the zero-temperature case), acting on 
the thermal state of the dual Hamiltonian. 

However, it does not appear possible to construct a 
Markovian effective noise model for nonzero temperature 
using the same techniques as for zero temperature. The 
reason is that our arguments were based on the assump- 
tion that the reduced state pk on each of the non-trivial 



sites in the dual state does not differ greatly from its 
value in the dual of the unperturbed resource state. This 
is indeed the case for small local perturbations to the 
Hamiltonian (as we prove in Appendix [D]) , but it need 
not be true for nonzero temperature. For example, con- 
sider the one-dimensional Ising model, with Hamiltonian 



i+l 



(71) 



(we include the Z 1 term to select out a unique ground 
state). In this model, it can be shown (e.g. using the 
transfer matrix method) that the reduced state on a sin- 
gle spin changes discontinuously as soon as the temper- 
ature is switched on (this is closely related to the disap- 
pearance of the magnetic order in the 1-D Ising model 
at nonzero temperature). Given the structure of the 
dual Hamiltonian as discussed in Appendix [B] there is 
good reason to believe that it will exhibit a similar phe- 
nomenon. 

The difficulty of treating thermal states in our frame- 
work should not be surprising, as the dual process has the 
perfect operation of the identity gate built in, whereas 
the cluster model is not expected to have a long-range 
identity gate at nonzero temperature. On the other 
hand, there exists a measurement protocol for a three- 
dimensional cluster model which retains the perfect op- 
eration of the identity gate at sufficiently small nonzero 
temperatures [41] . Therefore, if the dual process descrip- 
tion could be extended to measurement protocols such as 
this one, then it might be expected that the dual Hamil- 
tonian would possess an ordered phase that persists at 
nonzero temperature, such that the local reduced state 
varies continuously with temperature up to the phase 
transition. 



IV. TWO-DIMENSIONAL SYSTEMS AND 
FAULT TOLERANCE 

The equivalence we demonstrated in Sections [H] and 
TEH between MBQC on perturbed resource states and 
noisy quantum circuits, opens up the possibility of ex- 
ploiting the results in the literature on fault-tolerant 
quantum computation with noisy quantum circuits. 
Here, we will extend the results of the previous sections 
to the 2D cluster model, which, unlike the ID models 
considered previously, is a universal resource for quan- 
tum computation. We will again find that, provided the 
perturbation to the Hamiltonian respects a certain sym- 
metry, MBQC using the perturbed ground state as a re- 
source is equivalent to a noisy quantum circuit. We will 
show that the noise in this effective circuit description has 
no correlations in time (as in the previous section), nor 
any correlations in space. This reduction to local, Marko- 
vian noise will allow us to invoke the threshold theorem 
to deduce that, provided the perturbation respects the 
symmetry and is sufficiently small, the perturbed ground 
state remains a universal resource for MBQC. 
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1-D cluster state 



Uncoupled qubit 
(in |+) state) 



FIG. 7. The first step in generalizing the ID results to 2D 
models involves consideration of a 'quasi-ID' model, which 
consists of a ID cluster Hamiltonian acting on each of N 
qubit chains arranged in the vertical dimension, as well as 
a term favouring the |+) states on the uncoupled qubits. The 
model has a (Z2 x Z2) xN symmetry, arising from the Z2 x Z2 
symmetry associated with each of the N chains. We can treat 
this model as 'quasi-ID' by defining our sites (shown as green 
shaded areas) so that they span the vertical dimension. 



It should be emphasized that, although we make use 
of the theory of fault-tolerant quantum computation, our 
final result cannot be described as a fault tolerance result 
for MBQC, since it applies only to symmetry-respecting 
perturbations, and we assume noiseless operation of the 
measurement protocol. 



A. The 'quasi-ID' model 

Here, we make a first attempt at generalizing the ID 
results to a 2D model which is universal for quantum 
computation. The ground state of the model we intro- 
duce here is not strictly a universal resource for MBQC 
unless we allow non-single-qubit measurements; however 
the discussion here will serve as a stepping stone for con- 
sideration of the 2D cluster model in Sec. IIVBI 

In the absence of perturbations, the 2D model we con- 
sider involves N uncoupled ID cluster states arranged in 
the second dimension, as shown in Fig. [7J The Hamil- 
tonian acting on each chain is simply the ID cluster 
Hamiltonian. For generality we also assume the exis- 
tence of some uncoupled qubits, each with an associated 
term —X in the Hamiltonian (i.e. the ground state is 
|+)). In order to treat this 2D model within the same 
framework which we have developed for ID systems, we 
will consider an entire N v x 2 block (where N v is the 
extent in the vertical direction) to be a single 'site', as 
shown in Fig. [TJa); hence we can consider the lattice to 
comprise a ID chain of such 'sites'. The unperturbed 
ground state, which we denote I'Fq), then has an MPS 
representation which is essentially a tensor product of 
several copies of the ID cluster state MPS representa- 
tion, with a correlation system comprising N qubits. 




FIG. 8. In the 'quasi-ID' model, two-qubit gates in correla- 
tion space cannot be done with single-qubit measurements. 
However, after applying controlled-Z gates between neigh- 
bouring qubits in order to couple two chains, an entangling 
gate can be performed in correlation space by means of single- 
qubit measurements. 



Each chain contributes a separate Z2 x Z 2 symmetry, 
so that the model is invariant under a symmetry group 
G = (Z 2 xZ ? )* N = {($!,••■ ,g N )\gi,--- ,g N eZ 2 xZ 2 }. 
The projective representation of this symmetry in corre- 
lation space is the TV-qubit generalization of the Pauli 
representation, namely 



V((fli,--- >9n)) = V v {g{)®---®V v {g N ), 



(72) 



where Vp is the single-qubit Pauli representation of 
Z 2 x Z2, given by Eq. (10). It can be checked that this 



projective representation is maximally noncommutative, 
and because it has dimension 2 N = y/\G\, it must be the 
unique irreducible projective representation correspond- 
ing to its factor system (by Lemma [I] from Sec. II C I. 

Now, it is easy to see that, where S is the set of gates 
which can be executed in correlation space by a single-site 
measurement (up to Pauli byproducts) in the ID cluster 
state, we can execute any tensor product 



Si <K> s 2 



(g> sn, si, ■ ■ ■ , sn € S 



(73) 



(up to Pauli byproducts) in correlation space by a single- 
site measurement in our 2D model; we just do the cor- 
responding measurements on each chain separately. We 
can also find a measurement basis for a columnar site 
which induces entangling gates between two qubits in cor- 
relation space; however, this measurement basis clearly 
cannot correspond to single-qubit measurements, since 
the two chains would then remain uncoupled. For rea- 
sons that will become clear when we consider the relation 
of the present model to the 2-D cluster state in Sec. |IVB| 
we will only consider entangling gates between nearest- 
neighbour qubits in correlation space, for which we con- 
struct the measurement basis in a particular way, as fol- 
lows. 

We define the on-site unitary u, which involves ap- 
plying controlled-Z gates between neighbouring qubits to 
turn our original resource state \^q) into another graph 
state I'F'g) in which the two chains of interest are cou- 
pled, as shown in Figure [8] If A[-} is the MPS tensor 
for \^q) at the given site, then A'[-} = A[u^(-)] is the 
MPS tensor for \^'q). Using the measurement sequences 
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described in [32] , it can be shown that there exists a mea- 
surement basis {la)} for a columnar site, corresponding 
to single-qubit measurements, such that A' [a] = B a U, 
where U is an entangling two-qubit gate, and the B a 
are outcome-dependent Pauli byproducts. It follows that 
this two-qubit gate can be performed in correlation space 
(up to the same Pauli byproducts) by measuring in the 
basis {u^|a)}. 

From the above considerations, we see that the model 
we are discussing can be considered as a generalization 
of the 1-D cluster state in which N qubits can be propa- 
gated in correlation space, acted on by entangling gates 
between nearest neighbour qubits as well as single-qubit 
gates. In the presence of symmetry-respecting pertur- 
bations to the Hamiltonian, the arguments of Sections 
[TT| and |III| can still be applied for any finite N. How- 
ever, if we want to exploit the locality of the perturba- 
tion in the vertical direction as well as the horizontal, 
we need to make some additional arguments. First, we 
observe that (by Lemma [T]) the protected subsystem of 
correlation space (which corresponds to the ancilla sys- 
tem appearing in the dual picture of MBQC) will have 
dimension 2 , and by identifying the action of the sym- 
metry within the protected subsystem with Eq. (72 1, we 



can decompose the protected subsystem into N qubits, 
one associated with each chain. 

Our argument now hinges on two observations. First, 
the dual Hamiltonian of which the dual state is the 
gapped ground state, as constructed in Sec. |III C[ is in 
fact a sum of interactions that are local on the original 
two-dimensional lattice. Second, the unitary couplings 
Qk appearing in the dual process, which a priori could 
couple an entire columnar site to the entire iV-qubit an- 
cilla system, in fact acts trivially outside an appropriately 



localized area (see Fig. 10). These observations both fol- 



low from the form of the interaction 



z>i=5>><ii®n»i) 



(74) 



between a columnar site and the ancilla system. (Re- 
call that T>i and its inverse appeared in the development 
of the dual picture in Sec. |IID| as well as in the con- 
struction of the duality transformation T> from which the 
dual Hamiltonian H was obtained in Sec. IIIC ) It is 



easily seen that in the present quasi-lD setup, T)\ simply 
comprises a number of applications of the corresponding 
operator 1>i for the one-dimensional cluster chain (see 
Fig.|§. 

Now, consider a quantum circuit comprising a sequence 
of gates, and let Qk be the sets of physical (not ancilla) 
qubits acted on by the corresponding couplings Qk in the 
dual process. Just as in the one-dimensional case, we 
expect that if R = minfe 1> fc 2 dist(Qfc 1 , Qk 2 ) is much larger 
than the correlation length £ for the dual state, then the 

reduced state ^((j k Q k )" on [JkQk should be ap- 

proximately a product state over the Q^'s. Thus, arguing 
as in the one-dimensional case (see Sec. |II F I, we find that 



(1) 



(1) 



L 



FIG. 9. The unitary operator T> 1 , which couples a columnar 
site and the iV-qubit ancilla system. 



Qk 



(a) Single-qubit gate 



(b) Two-qubit gate 



FIG. 10. The unitary couplings Qk appearing in the dual 
process corresponding to single- and two-qubit gates. 



performing the measurement sequence on a perturbed re- 
source state corresponds to a noisy quantum circuit, with 
the noise described by a noise superoperator £ k following 
each gate. Furthermore, £k acts non-trivially only on the 
same qubits that were acted on by the corresponding gate 
in the original noiseless quantum circuit. The strength 
of the noise, as given by \\£k — ^27 1 1 <> , is determined by the 
deviation (in the trace norm) of the reduced density op- 
erator on Qk from its unperturbed value [see Eq. (49) in 
Sec. |II F) , which should be small for small perturbations. 

Let us now estimate the required scale-up in the size 
of the resource state. We only consider in detail the case 
of local quantum circuits (i.e. containing only gates act- 
ing between nearest-neighbour qubits). As in the one- 
dimensional case (Sec. HID), according to the rigorous 
factorization result proved in Appendix [C] the minimum 
separation i? m i n between any of the Qfc's required for an 
accuracy e scales like 



JW£ = O(m) + Opog(l/c)]. 



(75) 



The required scale-up can be expressed in terms of 
-Rmin, as follows. First, we must ensure that, at each time 
step, all non-trivial gates are separated by a distance of at 
least i? m i n . This leads to a scale-up by a factor of ~ i? m ; n 
in the number of time steps. Then, the buffering between 
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horizontal locations at which nontrivial gates take place 
implies another factor of R m i n scale-up. Hence, the total 
scale-up factor is s ~ R^in- ^ n ^ ne other hand, if the 
quantum circuit that we want to simulate is not already 
local, then translating it into a local circuit will introduce 
additional overhead (still scaling at worst polynomially 
in the number of qubits in the quantum circuit). 



B. The 2D cluster model 

Now we will return to the model we are actually inter- 
ested in: the 2D cluster model on a square lattice. In- 
vestigations of the effect on this model of perturbations 
[TTl 03H15] have demonstrated a variety of results de- 
pending on the perturbation. Here, we will focus on per- 
turbations respecting an appropriate symmetry. When 
this symmetry is enforced, the cluster model lies in a 
robust SPT phase, within which the identity gate is pro- 
tected and the effective noise model construction of this 
paper can be applied. 

In order to achieve our goal, we will establish an equiv- 
alence between the 2D cluster model and a 'quasi-lD' 
model of the type considered in the previous section. The 
basic idea is to define a (local) duality transformation U 
(not the same as the duality transformation T> which we 
have considered previously) which relates the two mod- 
els. Specifically, we define 



U 



n (° z ^ 



(76) 



(i,j)eL 



where (CZ)ij is the controlled-Z gate acting on qubits 
i and j, and the product is over an appropriate set L 
comprising nearest-neighbour pairs of qubits. By an ap- 
propriate choice of L, we can ensure that applying U to 
the 2D cluster Hamiltonian turns it into a model of the 
type we considered in the previous section. 

Now, let H be a perturbation to the 2-D cluster Hamil- 
tonian. Then UHW is a perturbation to the quasi-lD 
model, and the arguments of the preceding section can 
be applied provided that the perturbation respects the 
appropriate symmetry. Furthermore, the result (in terms 
of statistics of measurement outcomes) of performing the 



adaptive measurement protocol described in Section IV A 
on the ground state of UHU\ involving measuring the 
observables 6, must be the same as the effect of per- 
forming the same protocol on the ground state of H, but 
measuring the observables U^aU. We will now examine in 
detail this corresponding measurement protocol for per- 
turbations of the 2-D cluster model. 

In the quasi-lD resource state, there is a set of 'redun- 
dant' qubits which never need to be measured. It turns 
out to be convenient to assume, however, that we do 
measure those qubits, in the z basis, and that we do this 
before any other measurements. We observe that all the 
measurements that are performed on the quasi-lD state 
(after applying controlled-Z gates to couple chains where 



we want to perform a two-qubit gate) are all single-qubit, 
and are either in the z basis, i.e. measuring Z, or in the 
x-y plane, i.e. measuring tjg = (cos9)X + (sin9)Y for 
some angle 8. Hence, the corresponding observables to 
measure in the 2D cluster state are either of the form Zj 
(for some qubit j), or 



n ^ 



(77) 



for some qubit j and angle 6j , and where TV} is some set of 
neighbouring redundant qubits. But, since we measured 
the redundant qubits first, they are all now in eigenstates 
of Z. Therefore, labelling the measured values of Z on the 
redundant qubits by {zj}, we see that measuring Udjli^ is 

equivalent to measuring Og^ (^Y\keAf Zfc ) > wmcn m turn 

( j ) 

is equivalent to measuring a g and reinterpreting the 
measurement outcomes based on the value of IlfceAT z k ■ 
Therefore, we have shown that the measurement proto- 
col on the 2D cluster state can be implemented using 
only single-qubit measurements and adaptivity. It can 
be checked that the measurement protocol so constructed 
is essentially the same as the usual one for the 2D clus- 
ter state on a square lattice, which is described, e.g. in 

mm 

Finally, let us discuss the required symmetry. The du- 
ality transformation U can be used to relate the (Zi x 
Z 2 ) xN symmetry which protects the quasi-lD model to a 
corresponding one in the 2D cluster model. The form of 
the generators of the latter symmetry is shown in Figure 
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a) . Let us remark that we can also make similar argu- 
ments in the case that the ID chains are arranged on the 
2D square lattice in an unconventional way, for example 
diagonally [UJ [46] as shown in Figure [lljb). The ad- 
vantage of the diagonal layout is that the symmet ry [ one 
of the generators of which is depicted in Figure 11 [b)] 



takes a particularly simple form, due to the fact that ev- 
ery non-chain qubit neighbours an even number of chain 
qubits, and so the Z's that would normally appear on 
non-chain qubits all cancel. In particular, this symmetry 
commutes with an especially simple and physically mean- 
ingful perturbation, namely a uniform magnetic field in 
the x direction, i.e. V = B'^2 i X i . (The effect of such 
a perturbation has been studied numerically in [UJ HZ] ; 
the SPT cluster phase persists up to a first-order phase 
transition at \B\ = 1.) 



C. Perturbed ground states are universal resources 

Let us summarize the conclusions which we obtain 
from the considerations in Sections |IV A| and |IVB| by 
stating them as a theorem. In combination with the 
threshold theorem of fault-tolerant quantum computa- 
tion, this theorem will allow us to deduce that sufficiently 
small symmetry-respecting perturbations to the 2-D clus- 
ter Hamiltonian retain ground states which can serve as 
universal resources for MBQC. 
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FIG. 11. Two possible layouts for the ID chains (red lines) 
on a 2D lattice. Each layout is associated with a measure- 
ment protocol for MBQC on the 2D cluster state, and with a 
symmetry group (a representative generator of which is shown 
in each case). For any layout, we can construct an effective 
noise model corresponding to performing the associated mea- 
surement protocol on a perturbed cluster state, provided that 
the perturbation respects the associated symmetry. 



We consider perturbations to the exact cluster Hamil- 
tonian Hq on the 2-D square lattice, which we can take 
to be a sum of local commuting terms, with an energy 
gap to the first excited state of 2. Suppose now we con- 
sider a perturbed Hamiltonian H = He + V, and V is a 
perturbation of the form 



V 



V u , 



(78) 



where A is the set of all lattice sites, and each V u is an in- 
teraction term supported on the set B(u, r) of sites within 
some fixed distance r (more generally, interactions de- 
caying exponentially with distance would not present an 
obstacle to our arguments). We define the local strength 
of the perturbation by 



J = max II VU 



(79) 



The cluster Hamiltonian He belongs to a class of Hamil- 
tonians for which it has been shown [48J that the gap is 
stable to local perturbations, i.e. there exists a threshold 
1] > (depending only on r), such that the gap of the 
perturbed Hamiltonian is at least 1, provided that J < r\. 

Let us assume that the perturbation V respects an 
appropriate symmetry group, constructed according to 



of the symmetry groups depicted in Fig. 111. Suppose 
then we want to use the perturbed ground state |$) to 
simulate a local quantum circuit containing N qubits, T 
time steps, and m gates, with the gates drawn from the 
gate set S comprising single-qubit rotations, a two-qubit 
entangling gate (as constructed in Sec. IV A), and the 



non-unitary RESTART gate (which corresponds to the 
reinitialization of a qubit). We obtain the result 

Theorem 4. Provided J < r\, we can find an appropriate 
measurement protocol on the ground state |\&) such that 
the final reduced state on the output qubits is e-close in 
the trace norm to the outcome of the quantum circuit, 
with added noise. In each time step t of the equivalent 
circuit process, the appropriate gates are applied, followed 
by a noise process described by a superoperator £ t . This 
superoperator can be written as a tensor product £ t — 
®A^t,A, where the product is over 'locations', i.e. sets 
of qubits coupled by a gate in the time step t ( each qubit 
not coupled by a gate in the time step t also counts as a 
location, but £ t .A — 1 in that case). Thus, the noise has 
no correlations in space ( other than those due to gates 
acting between qubits) or time. Furthermore, the noise 
operator £t,A oi each location and time is close to the 
identity superoperator in the diamond norm: 



\\St,A -Z\\o <cJ, 



(80) 



for some constant c (dependent only on r). The number 
of qubits measured n satisfies 



n = NT{(D{m) + 0[log(l/e)]} . 



(81) 



Proof. The only aspect that we have not previously dis- 
cussed is the bound Eq. (80 1. Following the same ar- 



gument as in the one-dimensional case (Sec. II F), we 
find [using the analogue of Eq. (49 1] that the deviation 
\\£t,A — 2T|lo is bounded above by~A x = \\p x ~ px,o\\i, 
where X is the set of qubits in the 2-D lattice that affect 
the operation of the gate in question, and px and px,o are 
the reduced states on X of the perturbed dual state |^) 
and the unperturbed dual state j^o) respectively. Phys- 
ically, it should be clear that Ax will be small for small 
perturbations; in Appendix [D] we demonstrate that, so 
long as J < rj, the inequality Ax < cJ holds for some 
constant c depending only on r. □ 

Now that we have shown that perturbations in the 
Hamiltonian correspond to noisy quantum circuits, we 
can invoke the threshold theorem of fault-tolerant quan- 
tum computation [18, 49, 50 . For our purposes, the most 
suitable version is Theorem 13 of Ref. [18], which we can 
state as follows: 

Theorem 5. Let us assume a noise model as described 
in Theorem^ Then there exists a threshold rj' > and a 
constant a such that, so long as \\£t,A — ^ZT| | < rj' for all 
A,t, then the following propeties hold. For any e > 0, and 
any local quantum circuit C made from gates drawn from 
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S ( with N qubits, T time steps, and to gates ), there exists 
another local circuit C with gates drawn from S, such 
that C with noise produces the same result (in terms of 
the probability distribution for the final readout, and up to 
an error e) as C without noise. The scale-up factors for 
the number of qubits, the number of time steps, and the 
number of gates are all bounded by (const.) x log"(m/e). 

Combining Theorems [4] and [5j we obtain: 



Theorem 6. Consider the perturbed model H = 
He + V as described above. Then there exists some 
threshold rj" = min{TO rf/c} > (depending only 
on r) with the following property. Provided that 
J < rj", then for any local quantum circuit C (with 
N qubits, T time steps, and to gates), with gates 
drawn from S, we can find an appropriate measure- 
ment protocol on the perturbed ground state 
such that the result is equivalent (in terms of the 
probability distribution for the final readout, and up 
to an error e which can be made arbitrarily small) 
to the outcome of the original quantum circuit. As 
to — \ 00 with e held fixed, the number of measured 
qubits n satisfies 



n < NT x O (m log 3 " to) . 



(82) 



This is sufficient to show that the perturbed ground 
states remain universal resources, which is Theorem [3] as 
stated in Sec. PA" it is the main result of this paper. 



V. CONCLUSION 

In this paper, we have developed a framework to char- 
acterize the effectiveness of measurement protocols for 
MBQC with SPT-ordcred ground states of quantum spin 
systems. This has allowed us to prove the universality for 
MBQC of the ground states of perturbed versions of the 
2-D cluster Hamiltonian, provided that the perturbation 
is sufficiently small and respects an appropriate symme- 
try. 

The type of SPT order that we have presented here 
is that present in one-dimensional systems, which is re- 
lated to a nontrivial factor system (also known as a 2- 
cocycle). It is for this reason that, in order to establish 
universality in two-dimensional systems, we had to treat 
them as 'quasi-one-dimensional' and assume an extensive 
symmetry group (Z2 x Z?) xN , which grows with the ver- 
tical extent of the system. For standard, non-extensive 
symmetries in two dimensions, SPT orders can be re- 
lated to 3-cocycles [5TJ [52], but it remains to be seen 
whether similar connections can be drawn between such 
two-dimensional SPT order and MBQC. 

Finally, we note that if MBQC in ground states of 
quantum spin systems is to be a robust form of quantum 



computation, then it must be possible in the presence 
of arbitrary (not necessarily symmetry-respecting) local 
perturbations to the Hamiltonian, as well as at nonzero 
temperature. Non-symmetry-respecting perturbations 
break the symmetry that is essential to our argument; the 
difficulty of extending our treatment to nonzero temper- 
ature was discussed in Sec. IIII El Nor have we considered 
the effect of non-ideal measurements, or of decoherence of 
the resource state taking place during the course of the 
measurement protocol. Therefore, it remains an open 
question whether fault-tolerant MBQC is possible with 
such imperfections. 
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Appendix A: The dual FCS 

A pure finitely-correlated state (pFCS) [151 120] is 
the thermodynamic limit of the translationally invariant 
MPS generated by a fixed MPS tensor A. The nature 
of the correlations can be expressed through the transfer 
channel 



A{a) = A[k]aA[k] 



(Al) 



(here the sum is over some basis {\k)} for the site Hilbert 
space; it can be shown that this definition of A is inde- 
pendent of the choice of basis). In its canonical form, a 
pure FCS is further characterized by the following prop- 
erties: 

(a) A is unital, i.e. A(l) = I. 

(b) There exists a density operator A such that A^ (A) = 
A. 

(c) Defining a to be the largest magnitude eigenvalue of 
A^ other than the one corresponding to the eigenvec- 
tor A, we have that |o| is strictly less than 1. 

The correlation length is then defined by £ = —1/ log \ a\, 
and the eigenvalues of A correspond to the entanglement 
spectrum obtained from a cut in an infinite chain. 

Now we restrict ourselves to pFCS generated by ten- 
sors A satisfying the decomposition Eq. (16 1. We de- 



fine the CPTP superoperator V g according to V g (cr) — 
[V(g) <8> l]a[V{g) ® I] 1 ; it can be checked that V g is a lin- 
ear representation (in the space of superoperators) of the 
symmetry group G, and that it commutes with A' for all 
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g G G. Therefore, since A is the unique eigenvector of 
A^ with eigenvalue 1, it must satisfy V S (A) = x(ff)A f° r 
some scalars x(g)- The fact that V g is trace preserving 
ensures that x(ff) = 1- Hence we find that A commutes 
with V(g) <X> I; it follows by Schur's Lemma that A fac- 
torizes as A = SI ® A for some density operator A, and 
where SI = E/y[Gj is the maximally-mixed state on the 
protected subsystem (recall that the a/|G| is the dimen- 
sion of the protected subsystem). It follows that there 
is a -y/|G|-fold degeneracy in the entanglement spectrum 
throughout the SPT phase, generalizing the 2-fold degen- 
eracy in the (Z 2 X Z 2 )-protected Haldane phase [55] . 

We can define the transfer channel corresponding to 
the dual FCS (generated by ^4) according to 



A(p) = ^2A[k]pA[k]l 



(A2) 



Observe that A unital implies that A is also unital, and 
that 



A* (SI® a) = Sl®A\a) 



(A3) 



for any operator a acting on the junk subsystem. Hence, 
any eigenvalue of A must also be an eigenvalue of A. 
From this we can see that the dual FCS is also a pure 
FCS, and A is the unique fixed point of AJ; thus, the 
entanglement spectrum of the dual state is the same as 
that of the original state, but with the y / [G|-fold degen- 
eracy removecrl The respective correlation lengths obey 
the inequality £ < £. 

We now outline how the arguments of Section [II D| can 
be put on a rigorous footing within the pFCS formal- 
ism. Given everything that we have established so far, 
it can be shown that the reduced state p n of the origi- 
nal FCS on a block of n adjacent sites can be obtained 
from the corresponding reduced state p n of the dual FCS 
according to the quantum circuit shown in Fig. 
Assuming that we choose n large enough that al 
surements take place within this block of n sites, we can 
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then make arguments similar to those of Sec. II D and 
we find that the result of the adaptive measurement pro- 
tocol is equivalent to a sequence of interactions between 
the dual state and an ancilla particle as shown in Fig. 
12 b) , with the unitary interactions Qk defined as they 
were previously. 



Appendix B: Connection with the Kennedy- Tasaki 
transformation 

The Kennedy- Tasaki (KT) transformation [3U] is a 
non-local unitary transformation which transforms a 



2 A similar property was found numerically for the Kennedy- 
Tasaki transformation in |33| . We will discuss the connection 



between that transformation and our 
|Bl 



'dual state" in Appendix 



Pn 



(a) 



c> 



(b) 

FIG. 12. (a) The quantum circuit generating the reduced 
state on n sites for the original FCS. The ancilla particle is 
initially in the maximally-mixed state Q, and the symbol ■ 
indicates that it should be traced out at the end. (b) The 
adaptive measurement protocol acting on the original FCS 
is equivalent to a dual process acting on the dual state, com- 
prising a sequence of unitary interactions Qk coupling the dual 
state to an ancilla particle. 



spin-1 chain in the Z 2 x Z 2 symmetry-protected Hal- 
dane phase into a system where the symmetry is spon- 
taneously broken in the bulk. In this section, adapt- 
ing Ref. [33], we will define a generalized version of the 
Kennedy- Tasaki transformation, which can be applied to 
any system in the non-trivial SPT phase with respect to 
an on-site representation U(g) — [u(g)\® N of the symme- 
try group G = Z 2 x Z 2 (like the original KT transforma- 
tion, our generalization is defined for finite chains with 
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open boundary conditions and no terminating particles) . 
We will show that, when the ground state of the origi- 
nal system can be expressed as a pFCS, the ground state 
of the KT-transformed system is essentially the same as 
the state which we have referred to throughout this paper 
as the "dual state". We expect that for general ground 
states the situation should be qualitatively similar. 

Observe that the symmetry group is generated by two 
commuting elements x and z; hence for any g e G, 
we can write g = x m ^z n ^ for some m(g) and n(g) 
taking values of or 1. We will choose to write the 
unique non-trivial irreducible projective representation 
as V(g) = X m ^Z n( a\ where X and Z are the appro- 
priate Pauli operators. This is a rephasing of the Pauli 
representation Vp defined in Eq. (10 1; thus the factor 



system is different to, but in the same cohomology class 
as, that of Vp . (The construction of the dual state does 
depend on the specific choice of representative factor sys- 
tem for a cohomology class, although in a fairly trivial 
way; the present choice is the one that will ensure that 
the Kennedy- Tasaki transformation reproduces the dual 
state exactly.) 

The generalized Kennedy- Tasaki transformation T>kt 
is then defined as follows: 



V KT = Y[D kl , 



(Bl) 



k<l 



where D^i is a unitary coupling the two sites k and I 
according to 



D kl = \i){i\ ® u (a 



*(<») 



(B2) 



Notice t hat all the operators Dki appearing in the prod- 
uct Eq. (Bl ) commute. In the case that the particles are 



spin-l's, with the Z2 x Z2 rotation symmetry consisting of 
7r rotations about the x, y, and z particles (the simulta- 
neous eigenbasis of the symmetry for a single site is then 
{\x), \y), \z)}, where \a), a = x, y, z is the eigenstate 
of the spin component operator S a ), the transformation 
T>kt reduces to the standard Kennedy- Tasaki transfor- 



mation. Specifically, Eqs. (Bl) and (B2) correspond to 
Eqs. (4) and (5) in Ref. [33j. 

For open boundary conditions, we expect there to 
be an approximate (becoming exact in the limit as the 
length of the chain goes to infinity) four-fold degeneracy, 
due to the two-fold degeneracy associated with each edge. 
An appropriate analogue of the SPT pFCS for this choice 
of boundary conditions is as follows: the low-energy sub- 
space is spanned by states of the form 



\*(L,R)) 

IT- 



A 



A 



A 



A 



A 



R 



, (B3) 



where the MPS tensor A obeys the usual symmetry 
condition, where |L„) and are fixed end vectors, 

while \L) and (R\ are allowed to vary (within the two- 
dimensional irrep space) in order to generate the four- 
dimensional low-energy subspace. 

Now we want to examine what happens to a state of 
the form Eq. (B3) under T>kt- Suppose we set (R\ = 
(+|, \L) = |0). It is then straightforward to show [using 



the symmetry condition Eq. ( 12 ) and the decomposition 



Eq. ( 15 )] that applying all the pairwise interactions Du 
involving the first site gives the result 



JlAi) l*(o,+)> = 
\i>i / 



L* - A 



















+ 


A 




A 




A 




A 




R* 



(B4) 



Continuing this process, we find that 
2?jrr|*(0,+)) = 



.4 



.4 



.4 



.4 



.4 



, (B5) 



which is the dual state. To obtain the other states within 
the low-energy subspace, it is sufficient to observe that, 
using the symmetry condition on the tensor A, we have 



Mi,+)) = iu(x)r N \*(o,+)), 

\m-)) = M*)] 9N \m+))> 
Mi-)) = [ U (xz)]® N \m+)), 



(B6) 
(B7) 
(B8) 



and that T>kt commutes with u(g) for all g € Z2 x Z^. 
Hence, the other basis states for the Kennedy-Tasaki 
transformed system can be obtained from the dual state 
by application of a symmetry operation [u(g))® N for 
some g E G. The dual state thus represents one of the 
four distinct symmetry-breaking states within the low- 
energy subspace of the transformed system. 

Let us also note that the unitary transformation T> that 
we introduced in Sec. lIII 51 has much in common with the 
Kennedy-Tasaki transformation T>kt, even though they 
relate to different choices of boundary conditions. In- 
deed, an important property of T>kt is that for any local 
symmetry-respecting observable 6, T>ktoD kt remains 
local and symmetry-respecting. It turns out that this 
also holds for 2?, in the case of observables 6 acting in the 
bulk (the fact that T>oD> is local was established in Sec. 
Ill C it can be shown that T>oD^ still respects the on-site 
symmetry as well). Thus, like T>kt, the transformation 
T> can be applied to yield a local, symmetry-respecting 
Hamiltonian H in the bulk. This Hamiltonian need not 
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be the same as the Kennedy- Tasaki transformed Hamil- 
tonian, as can be seen by calculating how some simple 
symmetry-respecting interactions transform; specifically, 
for a chain of spin-l's, T> treats the x-, y- and z-directions 
on an equal footing, whereas T>kt does not. Neverthe- 
less, based on the results for pFCS described above, we 
expect H to have four degenerate, locally distinguish- 
able symmetry-breaking ground states in the bulk. On 
the other hand, when the transformation V is applied to 
the edge interactions (those which couple the ends of the 
chain to the terminating particles), the result need not 
respect the symmetry. Therefore, including the edge in- 
teractions favours one of the symmetry-breaking ground 
states over the others and leads to the non-degeneracy of 
the dual state |^). 



Appendix C: The factorization condition for ground 
states of a local Hamiltonian 



In this Appendix, we will show how to derive the 
approximate factorization result Eq. ( 69 ) for a non- 



The two-body result implies that the first term in the 



degenerate gapped ground state |\&) of a local Hamil- 
tonian. Our main tool is the existing theorem on the 
exponential decay of correlation functions for such a 
ground state [53"H55] . This theorem states that there 
exists a correlation length £ and a function f(x), with 
f(x) — C[exp(— as x — > oo, such that for any sets 
of lattice sites X and Y, and observables Ax and By 
supported on X and Y respectively, we have 

\(A x By)v - (Ax)v(B Y h\ 
\\Ax\\\\B Y \\ 

</(dist(A,r))min{|A|,|F|}, (CI) 

where denotes the expectation value of an observable 
with respect to |\E'), and \X\ denotes the number of sites 
contained in the set X. 

Suppose that instead of two sets of lattice sites, we 
have m disjoint sets X\,..., X m . Let N — Y^k=i l-^fel be 
the total number of lattice sites contained within all of 
the Xk's. We can obtain the following corollary: 

Lemma 3. For any observables Ax x , ■ ■ ■ , Ax m supported 
on the respective sets, 



\(A Xl ■ ■ ■ A x J* - {A Xl )y ■ ■ ■ (A Xm )v\ 



1-4 



\A 



x„ 



< f(R)N, 



where R is the smallest distance between any two of the 
X k 's, i.e. R ee min fe ^ dist(A fc , X t ). 

Proof. Without loss of generality, we can assume that 
\\A Xk \\ = 1 for all k. Observing that \(A Xl ) \ < \\A Xl \\ = 
1, we have that 

\(A Xl ■ ■•Axjv - (A Xl )* ■ ■ ■ (A Xm )*\ 

< \(A Xl ■ ■■ A x J* - (A Xl )^(Ax 2 '••Axjvl 
+ \(Ax 2 ■ ■ ■ A X J* - (A X2 )v ■ ■ ■ (A X JM- (C3) 



right-hand side of Eq. (C3) is bounded by \Xi\f(R). Con- 



tinuing in this way, the lemma follows by induction. □ 

Now we want to show that the reduced state on the 
union of the Xf-'s is close to a product state. To do 
this, we make use of the following lemma. We re- 
mind the reader that he re we make use of both the 
trace norm ||A||i = Tr V ATA and the spectral norm 
||A|| = max^|^,) = i ||A|^)||. They are both distinct from 
the norm induced by the Hilbert-Schmidt inner product. 

Lemma 4. Consider a vector space of dimension D. 



Then there exists a basis {Ei 



1,...,£> 2 } for the 



space of linear operators supported on the site, orthonor- 
mal with respect to the Hilbert-Schmidt inner product 
(A, B) = Tr(A^_B) 7 and comprising Hermitian operators 
Ei such that such that = 1. 

Proof. Given a basis {|ro)}, m = 1, • • • ,d, such a set 
can be constructed comprising the operators \m)(m\, and 
(l/Vf)(\m){n\ + \ri){m\) and (i/^){\m){n\-\n){m\) for 
m n. □ 

Now we are ready to prove the main result. 

Theorem 7. Let I'l') be the non- degenerate gapped 
ground state of a local Hamiltonian. Let p be the reduced 
state of 1^) on UfeLi -^-k, o,nd let be the reduced state 
on Xk ■ Then 



\\ P -p P ro d \\ 1 <Nf(R)d 2N , 

where p prod = (& k Pk- 



(C4) 



Proof. Recall that for a linear operator P, \\P\\i — 
maxo|| = i |Tr(AP)|, where the maximization is over all 
linear operators A with unit spectral norm. Now, we can 
expand 



aii > 



Ei, 



(C5) 



where the tensor product is over the sets X\, ■ ■ ■ , X m ; the 
the Ei are as constructed in Lemma [4] and the scalars 
oti lt ... t i m are determined by ctj^... j m = TrfA^iJ^ ® • • • <g> 
Ei m )], which implies (since ||A|| = 1) that 



Now, notice that 



< ll^lli-'-ll^lli- 



(C6) 



Tr(Ap)= J2 a ii,~,i m ( E h®-~® E iJp> ( C7 ) 



Tr(Ap prod ) = 



t (E h ) Pl ---(E im ) Pm . (C8) 
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Hence, by Lemma [3j we find that 
|TrL4(p-p prod )]| 

<Nf(R) ]T |ai 1 ,...,i ra |||S il ||...||£? i , 



(C9) 
(CIO) 



<Nf(R) (Il^ 1 lli'--Il^ m lli)(ll^ 1 l|---Il^ m ll) 



= Nf(R)d 2N , 
since ||Pfc||i||i?fc| = 1 for all k. 



(Cll) 
(C12) 

□ 



Appendix D: Local perturbations perturb 
continuously 

Physically, it should be clear that small perturbations 
in a gapped local Hamiltonian lead to small variations in 
the reduced state obtained from the ground state on a 
finite region of the lattice. Here we will give a rigorous 
proof of this fact, as follows: 

Theorem 8. Let H(s) be a differentiable path of Hamil- 
tonians of the form 



(Dl) 



ueA 



where the sum is over all the lattice sites u in a finite- 
dimensional lattice A, and H u (s) is supported on the set 
B(u,r) of sites within some fixed distance r of u. Sup- 
pose that for < s < S, the Hamiltonian H(s) has a 
unique ground state \^(s)), and there is a uniform lower 
bound 7 > on the gap. Then there exists a constant c 
(dependent only on the lattice geometry, on r, and on j) 
such that for any set X of lattice sites, we have 



\p x (S)-p x (0)\\ 1 <c\X\J / S, 



(D2) 



where p X {s) is the reduced state on X, i.e. p X (s) — 
Tr X c |#(s))(#(s)|, and J 1 = ma,x u eA, s e[o,S] \\d s H u (s)\\. 

Proof. The proof relies on the following consequence of 
the theory of quasiadiabatic continuation |56l I57j : under 
the given assumptions, there exists a family of Hamilto- 
nians T-L(s) such that 



*— 1*(*)> = w(*)|*(*)>, 

where "H(s) can be written as 

W( S ) = ^„( S ), 



(D3) 



(D4) 



such that for any site u, T-L u {s) can be approximated by 
an observable 'H u {s) supported on X c (the complement 
of X), with error 



(D5) 



where / is a rapidly decaying function (dependent only 
on r). 

Now, for s £ [0,S], we have (where p(s) = 

l*0OX*0OI) 



d a px(s)=iTr X c[H{s),p{s)}. 



(D6) 



Hence 



< 



< 



max \Tr(A x Tr X c[H(s),p(s)})\ 

Ax 


(D7) 


max \Tr (A x [H{s),p{s)])\ 

Ax\\=l 


(D8) 


max \Tr ([A x ,H(s)]p{s))\ 

Ax || = 1 


(D9) 


max \\[A x ,H(s)}\\. 

-4x||=l 


(D10) 


max y^\\[A x ,n u {s)}\\ 

Ax\\ = l 


(Dll) 



Here the maximization is over all linear operators A x 
supported on X with unit spectral norm. We have made 
use of the fact that for any linear operator P, ||P||i = 
max|| A || =1 |Tr(AP)|. 

Now, for any operator A x supported on X, we have 
that [A X ,T-L u (s)] = since A x and % u {s) are supported 
on disjoint subsets. Hence, using Eq. (D5) (and the fact 
that || = 1), we find that 



\[A x ,n u (s)}\\ < 2\\A x \\\\H u {s) - Hu(s) 

< 2 J7 ( dist( 7 "' X) 



so that 

\\d s p x {s)\\ 1 <2J'Y J f 



ueA 



min^gx dist(it, x) 

1 



We can bound the sum according to 
^ / dist( U ,AT 



ueA 



tiGA 



7 



m.m xeX dist(u, x) 



<^2^2 f ^ dist ( M ' x ) 



tiSA x£X 



EE/ 



/ dist(w, x) 



x<£X ueA 



< |X|max V / 



ueA 



7 

dist(u, x) 

7 



< \X\c/2, 



(D12) 
(D13) 

(D14) 

(D15) 

(D16) 

(D17) 

(D18) 

(D19) 
(D20) 



where in the last step the rapid decay of / ensures that 
the sum is bounded by a constant c/2 dependent on the 
lattice geometry, 7, and r. 
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Hence, by the triangle inequality for the trace norm, as required. □ 
we have 

\\px(S)-px(0)\\i< [ S Wd.pxie^da (D21) 
Jo 

< c\X\J'S, (D22) 
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